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The neutrinoless double beta decay (0u/3[3) is very interesting for several reasons. In 
the first place, this decay mode is viable only when the neutrino is a massive Majorana 
particle. As such, it constitutes a critical touchstone for various gauge models that go 
beyond the standard SU(2)l x £7(1) gauge model of electoweak interactions. Secondly, the 
neutrinos with nonzero masses have many interesting consequences for the history of the 
early universe, in the evolution of stellar objects, and the supernovae astrophysics. Thirdly, 
besides the issue of m v ^ 0, there are other open questions in neutrino physics the answers 
to which depend on Qvf3f3 decay, such as: Why does nature favor only left-handed currents? 
Does the majoron exist? Yet, we shall not understand the Ovftfi decay unless we understand 
the two neutrino double beta decay (2^/3/3) . The last one is the rarest process observed so 
far in nature and offers a unique opportunity for testing the nuclear physics techniques for 
half-lives ^ 10 20 years. Thus, the comprehension of the (3(3 transition mechanism cannot 
but help advance knowledge of physics in general. 

In recent years the quasiparticle random phase approximation (QRPA) has been the 
most popular method to deal with the problem of + — > + double beta decay [JTJ . Within 
this model the /3/3-decay amplitudes are extremely sensitive to the interaction parameter 
in the particle-particle (PP) channel, usually denoted by g pp . Independently of the nucleus 
that decays, of the residual interaction that is used, and of the configuration space that is 
employed, all the QRPA calculations done so far exhibit the following general features. 

(i) Close to the "natural" value for g pp (g pp = 1) the 2v/3/3 moments have first a zero and 
latter on a pole at which the QRPA collapses. 

(ii) The zeros and poles of the 0^/3/3 moments for the virtual states with spin and parity 
J 71 = 1 + are strongly correlated with the zeros and poles of the 2vf3f3 moments. 

(iii) The Ou/3/3 moments of multipolarity J 71 " ^ + , 1 + also possess zeros and poles but at 
significantly larger values of g pp . 

(iv) As a function of g pp , both the 2vf3f3 and Qv/3/3 moments always present similar 
shapes. 

Fig. [I] illustrates the behaviour of the + — > + (3(3 matrix elements for several nuclei. 
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In the upper panel the 2v(3(3 moments {M.2v) are shown. The other two panels contain the 
Ov/3/3 moments of multipolarity J 71 = 1 + (M. 0l/ (J n = 1 + )) and total 0v(3[3 moments (Aio v ), 
induced by the neutrino mass mechanism. These results have been obtained with a 5 force, 
using standard parametrization presented elsewhere ||10|| . Instead of the parameter g pp , I 
use here the ratio between the triplet and singlet coupling strengths in the PP channel, i.e., 
t = Vt/v s . Calculations with finite range interactions yield similar results f|-||. 

More that once we have pointed out that the (3(3 amplitudes go to zero within the 
QRPA because of the restoration of both the isospin and SU(4) symmetries. We have also 
suggested a physical criterion for fixing the PP coupling strength based on the maximal 
restoration of the SU(4) symmetry (t = t sym ). Yet, the general characteristics mentioned 
above suggest the existence of some additional regularities, and the present concern reflects 
upon a global understanding of the (3(3 transition mechanism within the QRPA. Only in 
this way one can get a full control of the calculations, which is one of the prerequisites for 
a reliable estimate of the nuclear matrix elements. 

To begin with, I resort to the single mode model (SMM) description |J of the /3/?-decays 
in the 48 Ca — > 48 Ti and 100 'Mo — > l00 Ru systems. This is the simplest version of the QRPA, 
in which there is only one intermediate state for each J w . 

In the SMM the Ou and 2v moments for the + — > + transitions read || 



where M\ v and Ai® u (J + ) are the corresponding unperturbed matrix elements. Here 
G(J + ) = G(pn,pn; J + ) are the PP matrix elements, uj° is the unperturbed energy, and 
ujj+ are the perturbed energies. I will assume that the isospin symmetry is strictly con- 
served, in which case Aio u (0 + ) = 0. This statement is also valid for full calculations and 
therefore no further reference will be made to the intermediate states J w = + . When the 
pairing factors are estimated in the usual manner, one gets 




(1) 
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uo = u°yfl + F(34 + 9F/u°)/25u° + 16G(1 + F /lo°)/25uj , (3) 



and 

u = + 4F(45 + F/uj°)/225uj + G(270 + 172F/o; + 49G/uj°)/225uj , (4) 

for the single pair configurations [Of 7/2(71)0 f 7 / 2 (p)]j+ in 48 Ca and [0g-//2(n)0gQ/ 2 (p)]j+ in 
100 Mo, respectively. Therefore, while the numerators in Eq. (|2]) depend only on the PP 
matrix elements, their denominators depend on the particle-hole (PH) matrix elements 
F(J + ) = F(pn,pn\ J + ), as well. The numbers in the last two equations arise from the 
pairing factors. As illustrated in Fig. the SMM is a fair first-order approximation for the 
2vl3f3 decays in 48 Ca and 100 Mo nuclei. 

The role played by the ground state correlations (GSC) in building up Eqs. @) and @) 
can be summarized as follows: 

(a) The numerator, i.e., the factor (1 + G/u°), comes from the interference between the 
forward and backward going contributions. These contribute coherently in the PP channel 
and totally out of phase in the PH channel. 

(b) The G 2 and F 2 terms in the denominator are very strongly quenched by the GSC, while 
the GF term is enhanced by the same effect. In particular, for 48 Ca the term quadratic in 
G does not contribute at all. 

It can be stated therefore that, within the SMM and because of the GSC, the 2v matrix 
element is mainly a bilinear function of G(l + ). Besides, it passes through zero at G(l + ) = 
—oj° and has a pole when uj\+ = 0. Similarly, all A4 0u (J + ) moments turn out to be quotients 
of a linear function of G(J + ) and the square root of another linear function of G(J + ). Both 
the zero and the pole of Aio u (l + ) matrix element coincide with those of the 2v moment. 
One also should bear in mind that the magnitudes of the interaction matrix elements G(J) 
and F(J) decrease fairly rapidly when J increases. Thus the quenching effect, induced by 
the PP interaction, mainly concerns the allowed Ov moment. For higher order multipoles it 
could be reasonable to expand the denominator in Eq. (H) in powers of G(J + )/u° and to 
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keep only the linear term. This term strongly cancels with a similar term in the numerator 
and the net result is a weak linear dependence of the M.q v {J + ^ 1 + ) moments on the PP 
strength. Obviously, for the last approximation to be valid, the parameter t (or g pp ) has to 
be small enough to keep U\+ real. Briefly, the SMM can account for all four points raised 
above, and leads to the following approximations for the dependence of the (3(3 amplitudes 
on the PP strength 

M 2 ^M 2u (t = 0) 1 -^/^, (5) 

L — t/ti 



and 



M u = Mou(r = l + ;t = 0) J /7/ " 



V 1 " t/h 

+ M 0u (r^l + ;t = 0)(l-t/t 2 ), (6) 

where t\ > t and t 2 ^> ti, and the condition t < t% is fulfilled. It is self evident that these 
formulae do not depend on the type of residual interaction, and that analogous expressions 
are obtained for the (3(3 matrix elements when the parameter g pp is used (with g pp, s for t's). 

The common behavior of the (3(3 moments for all nuclei, together with the similarity 
between the SMM and the full calculations for 48 Ca and 100 Mo (shown in Figs. [I] and 
respectively), suggests to go a step further and try to express the exact calculations within 
the framework of Eqs. (|5|) and (|6]). At a first glance this seems a difficult task, because: (i) 
the SMM does not include the effect of the spin-orbit splitting, which plays a very important 
role in the (3(3-dec&y through the dynamical breaking of the SU(4) symmetry, and (ii) the 
full calculations involve a rather large configuration space (of the order of 50 basis vectors). 
However, the reliability of formulae (|5[) and @ is surprising. The results are presented in 
Table |I[ In the upper, middle, and lower panels I show the values of the parameters to, tx, 
and t 2 that fit the (3(3 moments displayed in the same order in Fig. [I]. I also list the values of 
the moments A4 2u , M-q^P = 1 + ), and A4q u (J + ^ 1 + ) for t = 0, together with the quantity 



J\f = \JJ2t=o[-M-exact(t) — M.fi t {t)] 2 that is an index of the goodness of the fit. The largest 
error occurs for 100 Mo. Still, even here it is not possible to distinguish visually the exact 



curves from the fitted ones. (This makes needless the exhibition of the adjusted curves.) In 
fact, for this nucleus the proposed formulae reproduce better the exact (3(3 moments than 
those obtained from the SMM. It is also gratifying that all three fits yield quite similar 
values for to and t\. The differences are at most of the order of 10%. 

A comment regarding the full QRPA calculations might be appropriate. The matrix 
element M.2v can always be expressed by the ratio of two polynomials in C(l + ) and F(l + ) 
(see Eq. (8) of Ref. ||). For a n dimensional configuration space these polynomials are of 
degrees 2n-l and 2n, respectively. The above results seem to indicate that cancellations of 
the type (a) and (b) are likely to be operative to all orders, and that the linear terms in 
G(l + ) are again the dominant ones. General expressions for the 0u moments, as a function 
of the PP and PH matrix elements, are not known, but a similar cancellation may be taking 
place in these as well. 

In summary, I have designed the Eqs. (|5|) and (§) and verified that they nicely reproduce 
the full calculations of the (3(3 matrix elements evaluated with a zero range force. I also feel 
that they are of general validity, and that any modification to the nuclear hamiltonian or to 
the configuration space can only change the coefficients in these formulae, but will not lead 
to a different functional dependence. Thus, we possess now a global understanding of the 
(3(3 transition mechanism (and a full control of the calculations) within the QRPA, which 
was the aim of this letter. 

It should be stressed that for practical application one always has to perform the complete 
calculation in order to do the fit. The real advantages of the analytic formulas fl5|) and (^|) 
are: 

1) they exhibit, in a very simple way, the main physics of the /3/3-decay in the QRPA model, 
and summarize the common features of the calculations done until now, and 

2) they establish the potential and limits of the QRPA method, and give a hint of direction 
that should follow the future theoretical studies. 

The pole at t = t\ is the response of the QRPA to the nonphysical situation, in which 
the energy of the lowest virtual J 71 = 1 + state becomes = (Ei + Ef)/2, where E{ and Ef 
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are, respectively, the energies of the initial and final states. There is no reason in principle 
why this should not happen in a nuclear model calculation (for a sufficiently large value 
of t). But, within the QRPA approach the pole develops close to the "natural" value of 
t, which makes the (3(3 moments to vary rather abruptly in the physically relevant interval 
to > t ^ ti. Certainly, this is a weak point of the QRPA |TT| and it is not clear yet how it 
could be circumvented. 

A qualitative agreement, between the shell model and QRPA results for the 2vf3(3 matrix 
elements in A8 Ca, has been reported @,f|. When applied to medium and heavy nuclei, the 
shell model is always accompanied by a very severe truncation of the configuration space, in 
order to become tractable. Contrarily, the QRPA is a readily accessible and fully controlled 
approach, and as such it calls for further developments. Efforts in this direction have recently 
been done by extending the model to describe the 2v decays to an excited final state fL2| , 



and by including the core polarization corrections to the effective interaction [13 
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FIGURES 

FIG. 1. Calculated double beta decay matrix elements Miv (in units of [MeF] -1 ), 
Mou(J w = 1 + ) an d Mov, as a function of the particle-particle S = 1, T = coupling con- 
stant t. The 48 Ca nucleus has been evaluated within 2hu) and 3huj major oscillator shells. For 
the remaining systems I have adopted the oscillator shells 3hu> and Ahu plus the Ohg/2 an d O/17/2 
intruder orbitals from the 5hui shell. The "physical values" of the parameter t (t sym ) are shown in 
the last row of Table |. 

FIG. 2. The exact (solid lines) and SMM (dashed lines) matrix elements M.i v (in units of 
[MeU] -1 ), as a function of the coupling constant t/to (defined in the text). 
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TABLES 

TABLE I. The coefficients to, t\, and £2 and the matrix elements M.2u, Mou(J n = and 
■Mou(J n 7^ 1 + ) f° r t = 0, in the parametrization of the 2v and Of /3/3 moments. The quantity M 
is the norm of the residuals, i.e., the square root of the sum of squares of the residuals. The exact 
and fitted matrix elements are equal at t = 0, and the strength t is varied, by steps of At = 0.1, up 
to the collapse of the QRPA. The matrix elements M.2u are given in units of [MeE] -1 . The values 
of the PP coupling strength, which lead to maximal restoration of the SU(4) symmetry (t = t sym ), 



are shown in the last row. 
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